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^ I We show that the amphtude of the finite size correction term for the nth 

*^ ' moment of the partition function, for randomly interacting directed polymers, 

diverges (on the high temperature side) as {uc — n)~^, as a critical moment 
He is approached. The exponent r is independent of temperature but does 
depend on the effective dimensionality. There is no such divergence on the 
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The averages of therinodynainic quantities for quenched random systems are often done 
through the rephca trick which requires an analytic continuation of the integer moments 
of the partition function Z" to n — ;> (See e.g. |]I[), with the overbar representing the 
disorder average. It is only recently that several issues related to this continuation could be 
debated quantitatively because of the results available for the problem of a directed polymer 
in a random medium 0. This was not feasible for other models. For example, the replica 
prediction of the nonexistence of the higher order (> 3) cumulants of the free energy of 1+1 
dimensional directed polymer problem has been challenged by direct numerical computation 
of those 0,^]. In large dimensions, the behaviours of the moments for n > 1 and n < 1 are 
found to be different ^ . The possibility of the existence of a singularity in n has also been 
pointed out P,0|. All of these point finger to the analytic continuation and interchange of 
the n — > with the thermodynamic limit. In complex situations like the spin glass, even 
the effective hamiltonian for the nth moment can admit of many phases adding fuel to the 
firei. 

In this paper we show an even more striking possibility in a directed polymer system. It 
is a divergence of the finite size effect for the moments of the partition function in the high 
temperature phase as a critical moment is approached. This divergence is distinct from the 
known singularity, at the thermodynamic transition point, in the finite size correction to the 
free energy for disordered systems |jll| like spin glasses |p, the random energy model [10 



and directed polymer in a random medium ^. In contrast, the divergence we report here 
is not necessarily at the bulk transition point. 

Our model involves two directed polymers which interact on intersection at sites. The in- 
teraction energy is random with zero mean. The model is similar to those used for depinning 



transitions of line like objects [12]. In the continuum version, the random interaction model 



is known to exhibit a disorder induced transition from a weak disorder (high temperture) to 



a strong disorder (to be called the RANI phase) phase [T^JTj]- Fo^' detailed understanding 
of the phases, we study the model on a lattice using a real space renormalization group ap- 
proach. Since real space renormalization group can be implemented exactly on hierarchical 



lattices [jT5|, we consider the problem from the beginning on a hierarchical lattice. As usual, 
the procedure can be thought of as an approximate real space renormalization group on a 
regular lattice. 

The lattice, as mentioned, is generated hierarchically by replacing each bond at the /ith 
generation by a "diamond" of b branches (Fig 1). Two directed polymers start at the bottom 
of the lattice and meet at the upper end, without any backtracking. The model of interest is 
the site model where the two polymers interact whenever they meet at a site, with a random 
interaction energy e. For simplicity, all the directed paths are taken to be identical as far 
as the choice of e on these is concerned. However, the energies along a particular path are 
independent and uncorrelated, so that the randomness is only in the longitudinal (special) 
direction and not in the transverse directions (see fig 1). For this paper, the only information 
needed about randomness is the behavior of the moments of the weight y = exp(— e/T) (T 
being the temperature in units of the Boltzmann constant) which, for a gaussian distribution 
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of e, behaves as y™ = y^ . Other distributions have been considered, though will not be 
elaborated upon here. 

Let Zfj_ be the partition function for a given realization of randomness, and let S^ = h^f^^^ 
be the number of single chain configurations, at the /xth generation. Here L„ = 2^ is the 



length of a directed polymer. We define Z^{n) = Z"^/ SJ^, to factor out the free chain 
entropy, and call Z^{n) the moments. We establish that, like the random medium problem 
0, for a given temperature, there is a critical value n = nc{y) below which all moments are 
in their high temperature phase, in the thermodynamic limit of course. In this limit, Z^{n) 
approaches a fixed point value for n < ric, whereas, for n > ric, the moments diverge but 
with a finite "free energy" density ffj.{n) = {nL^)~^ In ^^(n) (RANI phase). The approach to 
the thermodynamic limit can be written generically as 

g,{n)=g{n) + B,{n)L-'^^^ + ..., (1) 

where g{n) is the thermodynamic limit (/i -^ oo), and Bg{n) is the amplitude of the finite 
size correction. We take 5'^(n) for n < ricto be [,Z^(n)]^/"' , and for n > rZc to be ^^{n). We 
will see that the corrections are power laws in L. 

One of the main results of the paper is the blowing up of Bz{n) as (n^ — n)~^, for 
n — > ric—, a feature whose existence was not even dreamed of before. No such divergence 
occurs for n > ric- We study the variation of r with the various parameters of the problem 
and establish its universality for a given distribution. We also characterize the behavior of 
f^(n) for n -^ nc+. 

For a given realization of disorder, the partition function can be written as (see Fig |1|) 

Z,^^ = bZl'^Wf^ + bib-l)Si. (2) 

The first term originates from the configurations that require the two directed polymers to 
meet at C, while the second term counts the "no encounter" cases. There are no energy 
costs at the two end points. The moments of the partition function, from Eq. ^, are 

Z^+,{n) = b-''J2Pnn.Zl{m), (3) 

m=0 

where Pnm = (m)(^ ~ i'^n~m:jjm:^ with the initial condition ZQ{n) = 1 for all the moments 
because there is no interaction in the zeroth generation (one single bond). 

The behaviour of the moments are determined by the stable fixed points of the recursion 
relation as /i ^ oo. The first moment has no fixed point for y > yi = 6^/ [4(6 — 1)]. For 
Yi = y/yi < 1, Z^{1) reaches a fixed point value Z*{1) for large /i. For yi > 1, f^(l) 
approaches a definite limit. For higher moments, the stable fixed points in extended spaces 
can be determined with a high precision on a computer by just solving a quadratic equation. 
It is easy to see that, for the stable fixed points of the first n — 1 moments, there is again a 
critical value ?/„, 

y-'=Ab-'-j:Pnm[Z*im)]\ 
m=0 

SO that for y„ = y'^/yn > 1 there is no fixed point for the nth moment eventhough the lower 
moments do have so. The stable fixed point, if it is real, can be written as 

Z*in)=b^{2r)-'[l-il-ynf^']. (4) 

The property to be used later is that 2Z*{n)y^b~^ < 1. 



Numerical analysis: By iterating the recursion relations for given b and y, we computed 
the moments of the partition function with large (70 digit) accuracy using Mathematica. 
Iterations upto 100 generations are done, and for checking, several cases with 300 generations 
are also considered. For a given y, there is a critical value ric so that for n < Uc the moments 
reach their fixed point value, as shown in the lower inset of Fig Q. This ric depends on y (i.e. 
temperature), b, and the distribution. For n < ric, the moments are in the high temperature 
phase, as per Eq. ^, because for these yn > y^- By the same token, for n > ric, the moments 
are in the low temperature phase. 

To analyse the finite size data (finite generations) for n < ric, we adopt the following 
numerical procedure. First, construct the differences A^(?7,) = Z^[^i{n) — Z}J^{n)^ from Eq. 
I, as 

A,(r^) = i?(n)(l-2-'^)L;^ (5) 

omitting, for simplicity, the subscript z. Hence, a log-log plot would give ifj and B{n), 
provided /i is sufficiently large. Such an analysis has been done for all the moments for 
various y and b. The sample plots of Fig ^ clearly show that the exponent ip is independent 
of n. It, however, depends on y. (The upper inset of Fig || shows a few sample plots 
for [-E^(n)]^/" — Z*(ny^'^ vs yU against B{n)L~'^ with the estimated values.) Fig ^ shows 
the growth of the amplitude with the moment indicating a divergence as the critical ric is 
approached from below. The location of n^ and the exponent r, as defined after Eq. ^ can 
be determined by choosing r such that B~^/'^{n) is a straightline with n. The intercept gives 

A similar analysis is done for f^(n) for n > ric. The exponent ■i/' = 1, and there is no 
divergence of the amplitude near ric- The finite size correction, from Eq. ^, is (Iny" — 
n\nb)L~^. In this case, the "free energy", f(n) in the thermodynamic limit vanishes in a 
singular fashion as Uc is approached . Assuming f(n) ~ {n — ncY, a and Uc can be estimated 
by linearizing the ^^l"{ri) vs n plot. 

The estimates of Uc from the two sides are consistent with each other. For n < ric, see Fig 
^, fits over a wide range in n give the exponent r as 0.71 ± .02 and this value is independent 
of temperature. We have checked this insensitivity to temperature upto Iny = 0.005 or 
Uc = 223.53. There is definitely a curvature near ric, and if only the last three or four points 
are taken, a lower value for r is obtained. We, however, believe that this or any other more 
elaborate procedure to estimate r is not warranted because the last few points are still far 
off from Uc- We, therefore, tend to believe the obtained value of r to be an upper limit like 
a mean field estimate. The free energy exponent a is found to be temperature dependent, 
i.e., nonuniversal. For example, for b = 4, a changes from 1.78 ± .03 at log^ = .04 to 
cr = 1.9 ± .05 at log^ = .08. What really happens at ric eludes us because it is not possible 
to hit an integer ric numerically. It is plausible that n — > nc±, and n = ric are to be treated 
separately. 

There is a strong dependence of the exponents on b and the distribution. For the gaussian 
distribution, the variation of r with b is shown in Fig |^. This dependence might be expected 
because as b is changed the effective dimensionality 26 = 2'^ changes. Interestingly enough, 
the exponent reaches a saturation as 6 — ^ 2, the minimum value of b for a transition. 

Linearized renormalization group analysis: To rationalize the results, we now analyze 
the linearized renormalization group transformation. 



For a given y, we linearize the recursion relation, Eq. ^ upto the moment for which the 
stable fixed point is reached (i.e., n < ric). Let us start with the situation where the moments 
are close to the fixed points, and define a column vector z of size no, the integer part of n^, 
with z{n) = 2^{n) — Z*{n) as the nth element. The transformation matrix R that takes the 
vector to a new one z' = Rz, is lower triangular with elements 26~"P„„i2*(ni), for m < n, 
and zero otherwise. R is of size no x no. 

The eigenvalues of this matrix R are just the diagonal elements A„ = 2Z* {n)y^b~"' , of 
which, in all cases, we find Ai to be the largest. By Eq. H, A^, < 1 for all n < ric- (See Fig |^ 
for the spectrum for a particular case.) If e„ is the nth eigenvector, then, after n iterations, 
2,M = ^^ AfcAj! e^, where A^ is the projection of the starting vector along e^. For large /i, 
it follows that the convergence to the fixed point for the nth moment is as BnL'"^, with 

^ = -log2 Ai = - log2 [l - (1 - yi)^/'] , (6) 

for all n, and amplitude Bn = AiSi „, where ei „ is the nth element of ei. 

There is a universality in the exponent ip when considered as a function of yi = y/yi. 
The numerical values for various temperatures, b, and distributions can be made to collapse 
on this curve. Incidentally, this is the same exponent that a pure system would have, with 
attractive interaction, if expressed as y/yc where yc is the binding-unbinding transition. 

The ratio Kn of n[2*(n)]*^"~^)/"i?„/ei „ is plotted against n in Fig ^ for the same situation. 
This shows the proportionality of the two. We have ensured that this proportionality is 
maintained for all the cases. Hence, we infer that the divergence of the finite size amplitude 
is really a consequence of the divergence of the elements of the eigenvector for the largest 
eigenvalue. 

How could the amplitude diverge? The eigenvalues are shown in Fig ^ for one temper- 
ature. In all cases studied, we find the rise of the last few eigenvalues approaching from 
below the first (largest) eigenvalue (Ai). This leads to an increase in the components of the 
first eigenvector e, whose ith component ei j is of the form 

where Ap = Ai — Ap. We, therefore, make the following hypothesis. In an analytic contin- 
uation in n, the largest eigenvalue of the transformation operator (not necessarily diagonal 
anymore) will be degenerate; the eigenvalues of the modes at the end, near n^, will come up 
and merge with Ai. It is this degeneracy that leads to the singular behaviour and produces 
a nontrivial exponent r. 

It is interesting to compare this situation with the problem of a directed polymer in a 
random medium [^^|l3 a^iid the similar randomly interacting directed polymers but now 
interacting on the bonds [1^. In these cases, unlike Eq. |], the temperature does not enter 
the recursion relations explicitly but only through the initial condition. For this reason, 
linearizing around the fixed point, with the associated eigenvalues and vectors, will not 
produce any strong n dependence, and hence no divergence in the finite size correction. In 
fact, an attempt (not serendipity) to formulate a problem, where the renormalization group 
transformations would contain all the important information, led us to this particular site 
version of the random interaction model. 



In the context of the random medium directed polymer problem, in 1+1 dimensions, a 
scaling form has been proposed for the moments, \ogZ^ = nLf + g{nL'^), where u is the free 
energy fluctuation exponent, and / is the thermodynamic free energy |^. In this case, the 
whole phase is the low temperature phase. In our case, for b > 2, there is a phase transition 
which is reflected through a nonzero Hc- The power law growth of the amplitude suggests 
a different scaling form for n < ric, namely logZ" = nLf + g{{nc — n)L^), where Q = ip/r. 
It is tempting to speculate that this fl is also the free energy fluctuation exponent. The 
question of scaling on the other side however remains an open question. 

To summarize, we have shown the existence of a diverging finite size effect with a weakly 
universal exponent, at a critical moment from the high temperature side - a scar left by the 
disorder. The thermodynamic limit on the high temperature side is identical to the pure 
system. The growth of the amplitude observed numerically can be a understood through a 
renormalization group argument. However, what controls these exponents remains a puz- 
zling issue. A theory to understand the universality (or its absence) of the exponent r is 
lacking. Full significance of the divergence is yet to be elucidated but it cautions that an 
analytic continuation in n has to be done with proper care. 

Partial support for this work has come from DST SP/S2/M-17/92. 
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FIGURES 

FIG. 1. Construction of a hierarchical lattice, (a) This is for b=2. Three generrations 
(;U=0,1,2) are shown. For this paper, identically marked points are taken to have the same inter- 
action energy e. But, e is random for points with different markers, (b) A more general motif with 
26 bonds. 

FIG. 2. Plot of log A^(n) vs logL^ for n = 5 (curve a), n = 20 (curve b), and n = 26 (curve 
c). These are for 6 = 4 and logy = 0.04. The straight lines are the fits to these log-log plots. 
The lower inset shows the plot of f(n) vs n. The transition is at Uc = 26.6. The lower moments, 
n < ric, are in the high temperature phase while the higher ones are in the low temperature RANI 
phase. The upper inset compares Eq. 1 with data points for n = 5(curve d) and n = 26 (curve e) 
on a log scale. 

FIG. 3. Plot of [B{n)]-^/'' vs n for 6 = 4 and (a) logy = 0.065, and r = 0.73, (b) logy = 0.04, 
and r = .73 (c) logy = 0.03, and r = 0.72 (d) logy = 0.02, and r = 0.72. The straightlines are the 
best fits through the points. 

FIG. 4. Dependence of the exponent r on 6. The line through the points is just a guide to the 
eye. 

FIG. 5. (a)The 26 eigenvalues of R for logy = 0.04 for which ric = 26.6. Note the log scale 
along y axis, (b) The horizontal straight line is the ratio Kn/Ki. This plot is on a linear scale, 
shown on the right. 



